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Controlling atomic and molecular processes by laser 
fields is one of current topics in physics and chemistry, 
and there are various control schemes applied to such 
processes. 1 These strategies are known to work when the 
system to be controlled is rather simple or small. How- 
ever, in reality, the system can become "complex," where 
the dynamics will be described by multi-level-multi-level 
transitions with a random interaction. Although the laser 
field can be obtained by optimal control theory (OCT) 1 
even for such complex systems, it is difficult to analyze 
the controlled dynamics because such a field obtained 
numerically is often too complicated to interpret. Hence 
it is desirable to have a more analytical point of view. 

It is well known that the it pulse 2 or its generaliza- 
tions can be employed to control few-level problems. 3 
Recently, an analytic result for multi-level control prob- 
lems between general quantum states has been reported. 4 
The scheme is based on STIRAP, 1 and assumes an in- 
termediate state coupled to the initial and target states. 
Though the scheme can accomplish perfect control, it re- 
lies on the energy level picture of a quantum system, so 
it is difficult to apply the scheme to large systems. 

In this short note, we propose a new approach to 
obtain an analytic optimal field for complex quantum 
systems. 5 Under a "coarse-grained picture" with OCT, 
which is valid for such a complex system, we derive an an- 
alytic expression for the optimal field which steers initial 
states to target states in a certain limit. By numerically 
solving Schrodinger equations, we confirm that perfect 
control is actually achieved. This point is important be- 
cause the zeroth-order solutions of OCT, 6,7 which look 
similar to our result, are not guaranteed to achieve per- 
fect control. Another point is that our final expression 
does not require a detailed information from the energy 
level picture, so that it is easy to apply to large quantum 
systems, in principle. 

We use OCT as a theoretical vehicle. The aim of OCT 
is to obtain an optimal field e(t) which guides the system 
from an initial state at t — to a target state |<I>/) at 
some specific time t — T. According to the OCT scheme 
by Zhu, Botina, and Rabitz, 1,8 the optimal field for the 
Hamiltonian H[e(t)] = Hq + e(t)V is given by 
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e(t) = —lm[m)\x(t))(x(t)\V\m] 
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where the quantum state \4>(t)) and the inversely- 
evolving quantum state |x(i)) satisfy usual Schrodinger 
equations with boundary conditions 1 0(0)) = and 
\x(T)} = |3>/)- In order to maximize the final overlap 
Jo = | (</ , (7 1 )| ( i > /) \ 2 , numerical iterations 8 are necessary. 

We approximately solve this problem by introducing 
"coarse-grained (CG)" Rabi states. We extend the usual 
Rabi oscillation between two eigenstates to the oscil- 
lation between two time- dependent states, \(f>o(t)) and 
\Xo(t)), defined by 

\Mt)) = U Q {t,0)\$i), \xa{t))=Uo{t,T)\$ f ) (2) 

where Uo(t2,t\) is the propagator for e(t) — from t\ to 
ti. We introduce the CG Rabi states as 

\4>(t)) = |0o(t))cosm-ie~ l9 |xoW)sin^, (3) 

\ X (t)) = -ie^oW) sin fi(t-T) 

+| X o(t)> cos Q(t-T). (4) 

The phase 9 is determined from a normalization condi- 
tion mmt)} = 1 [or ( X (t)\x(t)) = 1] as 

(Mt)\xo(t)) _ (MT)\*f) 
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(5) 



(xo(t)\Mt)) (*f\MT)) 

and f2 is a constant determined later. It is also shown 9 
that the CG Rabi states (3) and (4) satisfy Schrodinger 
equations with the optimal field derived below when 
and contain many eigenstates without a special cor- 
relation, i.e., they are random vectors, 5 and the target 
time T is long enough. 

Substituting eqs. (3) and (4) into eq. (1), and after 
some manipulations, we have an expression for the opti- 
mal field, 9: 10 

<t) = ^^Re [e^ s («M*)| V |xo W)] • (6) 

We have used K<^o(*)|xo(t))| -C 1 which is justified when 
and |$/) are random vectors as above. 
Substituting eq. (3) [or (4)] into the Schrodinger equa- 
tion with the field (6), and invoking the rotating-wave 
approximation (RWA), i.e., omitting highly oscillating 
terms 2 we have 



V 2 sin 2flT 



where 



V 




\(Mt)\v\ Xo (t)}\ dt 



(7) 



(8) 



is an "average transition element" between the time- 
dependent states (2). These equations constitute a set 
of equations to determine the value of f2 for given a and 
T. The final overlap follows as 

J = sin 2 QT, 



and the amplitude of the optimal field e is given by 




W)\ 2 dt = 
o v 



(9) 



(10) 
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Furthermore we set perfect control Jq = 1, which is 
realized when £1 = tt/2T [Sec eq. (9)]. This is equivalent 
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Fig. 1. The final overlap driven by the analytic optimal field, 
eq. (11), is shown as a function of the target time T. The system 
is random matrices with size N . The initial and target states are 
random vectors with Gaussian random elements, and the average 
Jo's for 100 ensembles are depicted. 
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Fig. 2. (a) The final overlap J = | (4>(T)\ */)| 2 and (b) the aver- 
aged field amplitude e are shown as a function of the target time 
T. Marks (x) represent numerical results by solving OCT equa- 
tions using the Zhu-Botina-Rabitz scheme. 8 Solid curves repre- 
sent our analytic results using the CG Rabi state. 



to take the limit a — > [See cq. (7)]. Substituting eq. (7) 
into eq. (6), and using ft = n/2T, we obtain the final 
expression, 

e(t) = ^Re[e- ie (Mt)\V\Xo(t))]- (") 
Note that this contains the same factor (<j>o(t)\V\xo(t)) 



as the zeroth-order solution of OCT, 6 ' 7 but the prefactor 
is introduced for perfect control. Actually the pulse area 
A(t) can be defined by 



A(t) 



/V) 
Jo 



2(Mt')\V\Xo(t')) 
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as an extension of the usual two-state case. 2 If we sub- 
stitute eq. (11) into eq. (12), and use eq. (8), we obtain 
A(T) = tt under RWA, i.e., eq. (11) is a generalized tt 
pulse. 

In Fig. 1, we show that the analytic field (11) actually 
works (Jo — > 1) for N xN random matrix systems. Since 
we have used the RWA and omitted 0(1/ y/N) terms 
in the derivation of cq. (11), perfect control is achieved 
when T, N — > oo as expected. In Figs. 2(a) and 2(b), we 
confirm eqs. (9) and (10) by comparing with numerical 
calculations using the Zhu-Botina-Rabitz scheme. 8 The 
agreement is good, especially for large T. 

We have studied random matrix systems as "complex" 
quantum systems, while real systems are located between 
simple systems and random matrix systems. Hence it will 
be interesting to apply this analytic field to more realistic 
cases of banded random matrix systems, i.e., a transition 
dipole matrix V is a sparse matrix and there can be no 
direct transitions between arbitrary quantum states. Our 
analytic field often fails to control in the situation that 
single-photon processes are forbidden, while we can show 
that nearly perfect control is achieved 10 if we restrict the 
initial and target state. Even in such a case, the next 
order solution is easily constructed along our picture, 10 
and can be applicable to multi-photon processes. 
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